Josephson current in diffusive multilayer superconductor/ferromagnet/superconductor junctions 



oo 

o 

O 

< 

in 

O 

o 
u 

Oh: 



O 

o 



ON 

m 
rn 

od 

o 

00 
O 



X 



Iver B. Sperstad/ Jacob Linder,' and Asle Sudb0' 

'Department of Physics, Norwegian University of Science and Technology, N-7491 Trondheim, Norway 

(Dated: Received August 25, 2008) 

We calculate the Josephson current in a diffusive superconductor/ferromagnet/superconductor junction, where 
the ferromagnetic region contains multiple layers (or domains). In particular, we study a configuration where 
there are two layers with an arbitrary relative in-plane magnetization orientation, and also include non-ideal 
interfaces and arbitrary spin-flip scattering. We study the O-tt oscillations of the critical current for varying 
junction width d, and find that the tt state vanishes entirely when the magnetic misorientation angle of the 
two layers exceeds a critical angle (j)c.- While cj>c tt/2 in the limit of high temperatures, we find that (fic 
becomes smaller than n/2 at low temperatures compared to Tc. O-tt oscillations are also found when varying 
the temperature or the misorientation angle for fixed values of d, and we present phase diagrams that show 
qualitatively the conditions for the appearance of such oscillations. We also point out how one may obtain 
significant enhancement of the critical current in such a system by switching the magnetization for selected 
values of the junction width d, and comment on the necessary conditions for establishing a long range triplet 
Josephson effect. 

PACS numbers: 74.20.Rp, 74.50.+r, 74.70.Kn 



I. INTRODUCTION 

The topic of superconductor-ferromagnet heterostructures 
has been a subject of intense research for several years^'^. 
Not only do such systems constitute model systems for in- 
vestigating the interplay between two fundamental condensed 
matter phenomena, ferromagnetism (F) and superconductivity 
(S), but recent advances in fabrication techniques of such hy- 
brid structures make applications increasingly attainable. Es- 
pecially S/F based Josephson technology holds great promise 
in nanoelectronics, e.g. as a physical realization of the qubit 
of quantum computation"^. Another possible device is in some 
sense analogeous to a spin valve exhibiting giant magnetore- 
sistance (GMR), i.e. strongly suppressing the cuiTent for op- 
posite orientation of the magnetization^ in two F layers sepa- 
rated by a normal metal(N). For our object of interest however, 
superconducting electrodes are used instead of ferromagnets 
and N is replaced with F, in which case the resistance effect of 
magnetization-switching is known to be reversed compared to 
the spin valval. 

The proximity effect between a superconductor and a nor- 
mal metal was predicted decades agc^ and has since been in- 
vestigated thoroughly both theoretically and experimentally. 
However, several new and interesting effects were predicted 
when the layer of normal metal was replaced with a ferro- 
magnet, due to spin-triplet correlations in the feiTomagnet 
induced by the exchange field^'''. Much attention has been 
given to SFS structures, which are studied as a somewhat 
more exotic class of Josephson junctions. The most interest- 
ing emerging phenomenon in SFS junctions is the appearance 
of the so-called 7r-state^'^, in which the difference in the su- 
perconducting phase across the junction is tt in the ground 
state, in contrast to the conventional state with phase differ- 
ence zero. The physical result of transitions between these 
states is usually a sign change in the critical Josephson cur- 
rent /c through the junction, the observable manifestation of 
which being a non-monotonic dependence of /c on parame- 
ters such as temperature and junction width. Been predicted 



for decades, the experimental verification of this phenomenon 
some years ago'"'" was one of the catalysts of the present ac- 
tivity on the field. Recently, the effect of magnetic impurities 
in SFS junctions was also investigated theoreticall y'^''^ . 

It is well known that in simple S/F structures, the proximity 
effect will only induce OSP triplet correlations (Opposite Spin 
Pairing, spin projection Sz = 0), and that ESP triplet cor- 
relations (Equal Spin Pairing, Sz — ±1) require inhomoge- 
neous magnetization^. ESP components are in some contexts 
referred to as LRTC components (Long Range Triplet Corre- 
lations), and are of special interest because they do not decay 
as rapidly in the ferromagnet as the other components, and 
may therefore evade the suppression of the supercurrent for in- 
creasing width of the ferromagnet. One way of achieving this 
in theory is to let the magnetization in the F layer have a he- 
lical structureJii^. A similar effect is also found by consider- 
ing a SIS junction where the S electrodes themselves exhibit a 
spiral magnetic order coexisting with the superconductivity'-. 
Yet another alternative is layered S/F structures'^ with non- 
collinear magnetization, or simply a SFS structure where F is 
replaced with several ferromagnets with different direction of 
magnetization. 

Recent theoretical studies have been focusing on junctions 
with two ferromagnets sandwiched between the supercon- 
ducting electrodes. There are two physical realizations of a 
system described by such a model. It may either describe 
a device with two distinct, consecutively placed F layers (a 
SF1F2S junction), constructed to achieve customized nonho- 
mogeneous magnetization. Several work a^- ' ^1 i^i^^i^ ^ i^^i^^ have 
considered the Josephson current in such heterostructures. On 
the other hand, the model may describe the more realistic ex- 
perimental situation of one ferromagnetic layer with several 
magnetic domains. Some models of the latter kind have in- 
cluded two or more in-plane magnetic domains2i24, and al- 
though this certainly is an interesting framework for under- 
standing real heterostructures, the present paper considers the 
two layers placed consecutively as a SF1F2S structure. The 
main motivation for this choice is that this configuration al- 
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lows for much easier experimental control of the magnetiza- 
tion. The misorientation angle may be tuned by applying an 
external weak magnetic field to the interface between the fer- 
romagnets if the magnetization axis is pinned in one layer, 
while in the other one there is an in-plane easy axis.- To ac- 
complish this experimentally, one would probably need some 
interlayer between the ferromagnetic films to avoid a locking 
between the corresponding magnetizations due to the interfilm 
exchange coupling. 

Calculations on the models referred to here have predicted 
O-TT transitions upon varying the strength of the magnetic ex- 
change field, the junction width or the temperature, depending 
on the relative orientation of the magnetization in the two F 
layer a^''^i'^'^'''^'i^^'^^ . For antiparallel magnetization, it is re- 
ported that the O-tt oscillations will vanish, rendering the tt 
state impossible^*. An enhancement of the critical current for 
the antiparallel orientation by increasing the exchange field 
was first reported by Bergeret et al.^, and shortly after elab- 
orated upon by others^^'^^. Much of the work has however 
been limited to the case of coUinear magnetization^^, but 
recently also SF1F2S systems with arbitrary misorientation 
angle for the magnetization have been analysed' '** '"^. In par- 
ticular, Crouzy et al. have shown' how the tt state of such 
a junction is suppressed for increasing misorientation angles, 
vanishing at a critical angle 0c ~ ti'/S when the temperature 
is close to the critical temperature, i.e. T/Tc ~ 1. One recent 
article-^ has even studied a corresponding ferromagnetic tri- 
layer structure, but focused chiefly on the LRTC contribution 
to the Josephson current in such a structure. It should be noted 
that while the majority of the work in this field is carried out 
in the dirty limit, considering diffusive F/S systems, several of 
the relevant paper a'^''^'^^'-^" study ballistic junctions as well. 

Ref. Jj points out the necessity of including additional ef- 
fects to get a more accurate description of such systems. The 
present article may thus be viewed as an extension of their 
work by studying a SF1F2S junction with noncollinear do- 
mains where non-ideal interfaces and magnetic impurities are 
also taken into consideration. For the latter, we will study the 
special cases where isotropic or uniaxial spin-flip scattering is 
present. Consequently, there are three questions addressed in 
this work which were not treated in Ref. Jj: i) how does spin- 
flip scattering influence the O-vr oscillations, ii) do nonideal 
interfaces change the qualitative behaviour of the system, and 
Hi) how does the Josephson current for such a system depend 
on the temperature. The possibility of investigating the lat- 
ter point is present in our model, as opposed to Ref. 1 , which 
was restricted to temperatures close to Tc- The reason for this 
is that the regime of weak proximity effect is only attainable 
either if the S/F transparency is low, or when transparency is 
high under the extra restriction that T « Tc, so that it is guar- 
anteed that the influence of superconductivity is weak in either 
case. 

For concreteness, we consider a diffusive Josephson junc- 
tion with two ferromagnetic layers with arbitrary in-plane rel- 
ative orientation of the magnetization, as shown in Fig. [T] Al- 
though we focus on this picture of distinct, controllable layers, 
the physically similar situation of magnetic domains will also 
be commented upon. The superconducting electrodes are two 



similar s-wave superconductors, and the interfaces between 
the superconductors and the ferromagnet are assumed to have 
low transparency. 
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FIG. 1: (Color online) The experimental setup proposed in this pa- 
per. Two s-wave superconductors are separated by two ferromag- 
netic layers with an arbitrary relative orientation of the magnetiza- 
tions. The ferromagnetic regions may model a domain-structure of 
a single ferromagnetic layer or correspond directly to two distinctly 
deposited magnetic layers. This is similar to the setup considered in 
Ref.m 

This paper is organized as follows. In Sec. |ll]we will briefly 
sketch the theoretical framework and go on to obtain a solu- 
tion for the proximity -induced anomalous Green's function in 
the ferromagnetic region of our SF1F2S system, from which 
an expression for the Josephson current can be calculated. In 
Sec. inilwe present the dependence of the Josephson current 
on the various parameters, analyse the resulting O-tt oscilla- 
tions in the junction and construct a corresponding phase di- 
agram. We discuss the applicability of our findings in Sec. 
IIVI and furthermore present a discussion of the absence of the 
long range Josephson effect in such SF1F2S systems. A sum- 
mary and some final remarks are given in Sec. |V] We will use 
boldface notation for 3-vectors, . ? . for 4x4 matrices, and 
for 2x2 matrices. 



11. THEORETICAL FORMULATION 

We address this problem by means of the Usadel equa- 
tion in the quasiclassical approximation. This corresponds to 
integrating out the dependence on the kinetic energy of the 
Gor'kov Green's function, obtaining thus the quasiclassical 
Green's function g as the object used to describe our sys- 
tem. This approximation is valid as long as all relevant en- 
ergy scales are much smaller than the Fermi energy ep, and 
correspondingly that all relevant length scales are much larger 
than the Fermi wavelength. The latter condition is reconciled 
with the presence of sharp interfaces in our model by intro- 
ducing appropriate boundary conditions, as discussed below. 
The approach is based on the Keldysh formalism for nonequi- 
librium superconductors, which is convenient to work in also 
in the present limiting case of equilibrium. Here one operates 
with 8x8 matrix Green's functions, in which the retarded 
Green's function , the advanced Green's function g^ and 



3 



the Keldysh Green's function are 4 x 4 matrix compo- 
nents. As both the advanced and the Keldysh component of 
the matrix can easily be expressed by the retarded component 
in our case, it will be implicitly assumed in the following that 
the Green's function under consideration is the retarded com- 
ponent . 

We will desist from further discussion of quasiclassical the- 
ory, and instead refer the reader to the considerable litera- 
ture that covers the Keldysh formalism and nonequilibrium 
Green's function a^''^^-^-^'-^^--^^ . We go on to write up the matrix 
structure of our quasiclassical Green's functions. In the bulk 
superconductors, the Green's function reads^ 



gs 



ic 



lT2Se 



-Ic 



(1) 



where c = cosh((?), s = sinh(0), and 6 = atanh(Ao/e), 
with Aq denoting the amplitude of the superconducting gap. 
The different signs of the phase x above correspond to the left 
(lower sign) and right (upper sign) superconducting bank, re- 
spectively. When not being in proximity to a superconductor, 
the Green's function for a bulk ferromagnet reads 



9F.0 = 



1 

-1 



(2) 



When being influenced by a superconductor, off-diagonal el- 
ements is introduced to this Green's function, and for weak 
proximity effect it is changed to gp « gp o + f. This pertur- 
bation can be expressed as^' 



/ 



lie' 
-[/(-£)]* 



(3) 



where the constituting anomalous Green's function can be 
written as a matrix in spin space on the forrn^^ 



/ 



/t /. + ./s 

ft - /s fl 



(4) 



Here, /s denotes the singlet component, ft the OSP triplet 
component and /| and fi the ESP triplet components, and it 
is these anomalous Green's function that the Usadel equation 
is to be solved for. 

In our calculations, we will account for the possibility of 
both uniaxial and isotropic spin-flip scattering by the parame- 
ter rjxy and rj^ as follows: 

Uniaxial spin-flip: 'qxy = 0, rj^ = 3, 

Isotropic spin-flip: ij^y = 1, ?72 = 1- (5) 

The spin relaxation time for spin-flip scattering will be de- 
noted Tsf, and is to be considered as a phenomenological pa- 
rameter in our approach. 

In the ferromagnetic regions Fi and F2, the linearized 
Usadel^^ equations take the form 

Ddlift ± /s) + 2i(e ± /i cos (/.)(/, ± /,) 
- ^ivzft ± 3/.) ±hsmcl>{f^ + /x) = 0, 
DdlU + (2ie - ^)fa - 2hsmcPf, = 0; a =t, i, (6) 



with (j) = Oin Fl. Eqs. ^ constitute a set of coupled, second- 
order, linear differential equations. Although an explicit an- 
alytical solution may be obtained for {/s, ft, f^} by solving 
Eq. (|6]l brute-force for non-zero (p, the resulting expressions 
are very large and cumbersome. We therefore proceed via an 
alternative but equivalent route. By a change of spin basis to a 
quantization axis which is aligned to the exchange field in F2, 
one obtains the equations 

DdUf: ± ./s) + 2i(£ ± h)if: ± /,) - ^{vj: ± m = o, 



Da|/; + (2i£-|j)/; = o. 



(7) 



The superscript ' denotes the new spin basis, and the s-wave 
component transforms as a scalar under spin rotations: — 
/s. The general analytical solution for these equations in the 
case of isotropic spin-flip scattering was obtained in Ref. [HI 
In the present case, we obtain 



I. 



-(ciK+e ''-^ + c2K^e'^-^ 



(8) 



where we have defined 



q± - V-DT,fi±p - 3 - r;, + 8ir,fe)/(2Z?r,f), 

K± = 3 - r]z ±p, p = y^(3 - r/z)2 - Gir^fh'^, 

k = ^[2ie-rixy/{2T,i)yD. (9) 

Taking 772 — rj^y ~ 1, the above expressions reduce to those 
of Ref. 38. Once the {/^ , /,', /^} have been obtained, one may 
transform them back to the original quantization axis || z. If 
we write 



/ = (/s + f -Z^Zi 



(10) 



one may from Eq. (01 identify the vector anomalous Green's 
function 



[/i-/T,-i(/T + /i),2/.]/2. 



(11) 



This is equivalent to the dk-vector formalism described in for 
instance Ref. [IB- For future use, we define f± — ft ± fs- 
Finally, the transformation to the new spin basis is (see also 
Fig.© 
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COS — sin ( 
, sin (b cos ( 



(12) 



In general, the linearization of the Usadel equation is a valid 
approximation in the case of a weak proximity effect. This 
may be obtained in two limiting cases: i) the barriers have 
low transparency or ii) the transparency is perfect (ideal in- 
terfaces) and the temperature in the superconducting reservoir 
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FIG. 2: The change of spin-basis from a quantization axis || z to a 
quantization axis || h. 



is close to Tc, such that Aq is small. An analytical approach 
is permissible in both scenarios, with differing boundary con- 
ditions. In case i), the standard Kupryianov-Lukichev (K-L) 
boundary conditions''^ are usually employed in the literature, 
while case ii) implies continuity of the Green's function and 
its derivative. In an experimental situation, the barrier region 
can hardly be considered as fully transparent, such that the K- 
L boundary conditions are more realistic than continuity of the 
Green's function and its derivative. We will therefore employ 
the K-L bondary conditions in this paper—. 

To obtain the anomalous Green's function, we must supple- 
ment the general solution in Eq. (O with the K-L boundary 
conditions at three interfaces. At the S/F^ interfaces located at 
X = —di and x = d2, one obtains 

27c^i5Fi9xffFiU=-di = [9s{-x),9Fi]\x=-di, 
'2jd29F2dx9F2\x=d2 = [9F2,9s{x)]\x=d2, 



We have defined the parameter 
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Rb 



(13) 



(14) 



where Rb is the resistance of the barrier region and Rp is 
the resistance of in the diffusive ferromagnetic regions (as- 
sumed to be the same for both Fi and F2). For the F1/F2 
interface, which denotes the separation of the ferromagnetic 
layers, we assume the resistance to be much smaller than at 
the S/Fi interfaces. Therefore, we model this by continuity of 
the Green's function and its derivative: 



9Fi = 5f2U=o, dx9Fi = dx9F2 



(15) 



Let us comment briefly on the case where Fi and F2 are two 
domains of a single layer A domain-wall resistance may quite 
generally be defined as = R — Rq, where R and Rq are 
the electrical resistances with and without a domain wall (i.e. 
homogeneous magnetization), respectively. When the width 
of the domain wall increases, i?^, 0, and vanishes all to- 
gether when the width of the domain wall is much larger than 
the Fermi wavelength"*'. In the present paper, we consider an 
abrupt change in magnetization at the interface of the two do- 
mains, corresponding to a very thin domain wall, such that 
one would in general expect a finite contribution to the resis- 
tance of the junction. To reduce the number of parameters in 
the problem, however, we assume that this resistance is much 



smaller than at the S/F interfaces and effectively set it to zero. 
In the case where Fi and F2 are separate ferromagnetic lay- 
ers, one may neglect the resistance at the interface between 
them by assuming a good electrical contact achieved during 
deposition of the layers. 

The {f2±; f2a} anomalous Green's functions are related 
to a a set of anomalous Green's functions in a rotated basis 
{/k, /LI via 

/2T = ^[cos0(/^^ + /^^) + 2isin(/./^^, + /^^ - /^J, 
f2i = ^[cos0(/^^ + /^^) + 2isin#^_, - /^^ + /^J, 



(16) 



where {/2-1-, /2ct} have the general form as shown in Eq. ([8]l. 
The complete anomalous Green's functions in the regions Fi 
may be written as 

/i± =feie-''--L±(0) + 62e«-"£^(0) 
+ 63e«+^L±(0) + M^'^+'^i^ (0), 



while in F2 one finds 

/2± = cie-'>--Lf{4,) + C2e'?--L±(0) 

+ '-^[iC,+C,)e^'^ + iD,+D 



(17) 



-ikx] 



f2a 



1. 



{cos0(CT+Cx)+a(CT-Cx)} 



2' 

+ 2isin0(cie-«-^ + C2e'?-^ + 036*+^ + C4e-«+^) 



(18) 

Above, we have defined (0) = cos (/)±'iK±/(8Tsf/i). Note 
that Eq. ( fTSl l reduces to exactly the same form as Eq. (fTTI i 
for (j) = (parallel magnetization), as demanded by consis- 
tency. The remaining task is to determine the 16 unknown co- 
efficients {bi}, {ci}, {Act, Ba, Ca, D„}. For clarity, we write 
out the boundary conditions explicitly. At a; — ~di, one has 

=T,i) 

ldidxfi± = c/i± T se""^, ■ydid^fia = cfia, (19) 
while at x — d2 we find 

^d2dxf2± = ±se'^ - c/2±, ld2dxf2a = -Cf2a, 



Finally, at x = one obtains 



/l± = /2±, dxfl± = dxf2±, 
fla ~ f2a, dxfla = dxf2a- 



(20) 



(21) 



Inserting Eqs. ([TtI i and ( fTSl l into Eqs. (fT9]l-(l2ni yields a set 
of linear equations which may be represented by a 16 x 16 
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matrix, and the solution for the 16 coeffiecients is found nu- 
merically. Once the anomalous Green's functions have been 
obtained, one may calculate physical quantities of interest. In 
the present paper, we will be concerned with the Josephson 
current 

/CO 
deRe{M+{e) + M-{e) 
-OO 

M|(e) + Mi(e)} x tanh(/3e/2), (22) 

with the definitions (cr =t, i) 

M,(£) = [fA-e)Vd^Ue) ~ /.(e)a, [/,(-£)]*, 
M±{s) = [f±{-e)rdj^{e) - /±(e)a,[/^(-e)]*. (23) 

The matrix ps in the first line of Eq. ( |22] | is defined by = 
ciiag(l, 1, —1, —1). The normalized current is defined as 

Iix)/Io = Miix,x)\/iNFeDAl), (24) 

which is independent of x for x G [— c?i, ^2] (due to conser- 
vation of electrical current). The maximal supercurrent the 
system can support, known as the critical current, is given by 
/(. = /( in the case of a sinusoidal current-phase relation. 

Before proceeding to disseminate our results, we briefly re- 
mind the reader (see e.g. Refs. 3.42) of the qualitative physics 
that distinguishes S/F proximity structures from S/N systems, 
and thus gives rise to e.g. O-tt oscillations of the critical cur- 
rent (which will be discussed in detail for our system in Sec. 
Ullt . The fundamental difference between the proximity ef- 
fect in an S/N structure as compared to an S/F structure is that 
the Cooper pair wavefunction acquires a finite center-of-mass 
momentum in the latter case due to the Zeeman-energy split- 
ting between the |- and |-spins constituting the Cooper pair 
The finite center-of-mass momentum of the Cooper pair im- 
pUes that the condensate wavefunction in the ferromagnetic 
region displays oscillations in space, permitting it to change 
sign upon penetrating deeper into the ferromagnetic region. 
Quite generally, one may write that the Cooper pair wavefunc- 
tion (order parameter) as = ^'oe"'^'"^ in an S/N structure 
while 'i! — 'fo cos(fc/.ia;)e~*''^'2a; jjj 5/p structure, where 
{kg, fc/.i, kf^2] are wavevectors related to the decay and os- 
cillation lengths of the proximity-induced condensate in the 
non-superconducting region. 

The fact that the proximity-induced superconducting order 
parameter oscillates in the ferromagnetic region suggests that 
the energetically most favorable (ground-state) phase differ- 
ence between the superconducting reservoirs might not al- 
ways be zero, as in an S/N/S junction. For a very thin fer- 
romagnetic layer, 'I' does not change much and there is no 
reason for why there should be an abrupt discontinuity in the 
phase at one of the F/S interfaces - hence, the system is in the 
0-state. If the thickness of the ferromagnetic layer is compa- 
rable to the oscillation length of 5" (~ in our notation 
above), then may cross zero in the middle of the ferromag- 
netic region and display antisymmetric behaviour. This is ac- 
companied with a shift of sign of the order parameter in the. 



say, right bulk superconductor as compared to the left bulk 
superconductor Under such circumstances, the energetically 
most favorable configuration corresponds to a phase differ- 
ence of TT between the superconductors, since A — Ao in one 
of the superconductors and A = — Ag — Aqb"^ in the other 
superconductor 

This is related to O-tt oscillations of the critical current as 
follows. The energy of the Josephson junction may in the 
tunneling limit be well approximated by 

ej-//,o(l-cos2x), (25) 

where /q contains the magnitude and sign of the critical cur- 
rent while 2x is the phase difference between the supercon- 
ductors. Now, // o is closely related to the proximity-induced 
condensate wavefunction 'J in the ferromagnetic region, and 
may change sign depending on for instance the width of the 
ferromagnetic layer Depending on its sign, the ground-state 
configuration corresponds to either 2% = or 2% = tt, and 
the critical current supported by the system will change sign 
depending on which of these phases the system is in (although 
the critical current itself is given by 2% = 7r/2). 



III. RESULTS 

Unless otherwise stated, we will fix /i/Aq = 10 and 
7 = 5 to model a realistic experimental setup with weak fer- 
romagnets (the exchange field was estimated to 5-10 meV in 
Ref. 43.) The particular choice of 7 is motivated by the fact 
that we expect the resistance of the barrier region to exceed the 
bulk resistance of the ferromagnets, and in addition to a low 
transparency of the S/F interfaces this clearly suggests that 
7^1. This will ensure a weak proximity effect, as explained 
in the previous section. In order to reduce the number of free 
parameters further, we in general set the widths of region Fi 
and F2 to be equal: di = d-z ^ d. The superconducting coher- 
ence length ^ — -y/D/Ao will be used as the unit in which the 
widths are measured. Where spin-flip scattering is included, 
we will use the parameter g = 1/ (rsf Ao) as a measure of the 
strength of this effect. 

We have numerically confirmed that the system follows a 
sinusoidal current-phase relation regardless of the direction of 
magnetization and all other variable parameters of our system. 
This allows us to focus on the state that supports the critical 
current, namely x = 7r/4. A sinusoidal current-phase rela- 
tion is moreover what should be expected for such systems 
with weak proximity effect"^-, and where the bulk supercon- 
ductors have the same symmetry We have also confirmed 
that the assumption of weak proximity effects holds by assur- 
ing that the value of the anomalous Green's function always 
obey \fs.t,a\ ^ 1 for the parameter range we consider. To 
make the computation of the solutions to the Usadel equa- 
tion numerically stable, we furthermore add a small imaginary 
term 16 to the excitation energy e, where the value S = lO^'^ 
has been used. This can be motivated as a way to account for 
inelastic scattering processes"*"*, interpreting the term as the 
inverse (positive) quasiparticle lifetime. 
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A. Zero temperature 

First, we will consider the case where the temperature is 
fixed to zero unless otherwise stated, i.e. the calculations are 
made with T/Tc = 0.001. The critical current as a function of 
the junction width d is shown in Fig. [3]for isotropic and uni- 
axial spin-flip scattering, respectively. Considering first the 
parallel case ip — 0, the well known O-vr-oscillations are re- 
produced, where the current change sign for certain values of 
d. It should however be noted that we have chosen to always 
plot the critical current as a positive quantity, as defined by 
Eq. (l24l l. because this is what is most commonly measured 
in experiments^''". We have confirmed that the oscillations 
are almost exactly periodic for the parameter range consid- 
ered here. Increasing the effect of spin-flip scattering tends to 
move the transition points between the state and the tt state 
towards higher values of d. Throughout our investigations, 
we find no significant difference between isotropic and uniax- 
ial spin-flip scattering for the width dependence of the critical 
current. Thus, we consider only uniaxial spin-flip scattering 
whenever the role of magnetic impurities is studied. 

The subplots of Fig. [3] show how increasing the relative 
angle of magnetization to (/) = tt/A shifts the first transi- 
tion points to the right and the second to the left, reducing 
the width region in which the system is in the tt state. For 
(j) — 7t/2, the oscillations have ceased entirely, leaving the 
junction in the state for all values of d. These effects are 
shown more clearly in Fig. HI which also shows that the oscil- 
lations do not return for (p e [7r/2, tt]. This can be expressed 
as a critical misorientation angle (f>c < 7r/2 over which the 
TT State is not realizable, which is in agreement with the find- 
ings of Ref. 1. Ref. 1 claims that (pc = 7r/2 independent of 
parameters as long as the system is near the critical tempera- 
ture, T/Tc ~ 1. We find (pc to have a somewhat lower value 
(j)c ~ 0.467r < tt/2 in the present case of T/Tc ^ 1, but we 
will show in Sec. IIIIBI how this value approaches tt/2 for in- 
creasing temperatures, and how it changes for other values of 
the exchange field than our particular choice of h/ Aq — 10. 
In light of the discussion concerning the qualitative physics 
involved in a S/F/S proximity structure, it seems reasonable 
to suggest that the vanishing of the O-vr transitions are directly 
linked to a strong modification of the Cooper-pair wavefunc- 
tion oscillation length inside the ferromagnet, which renders 
the vr-state inaccessible. Possible explanations will be dis- 
cussed further in Sec. |IV] 

As stated in the introduction, several works have contrasted 
the cases of parallel and antiparallel orientations, while the in- 
termediate angles have not been studied thoroughly (see, how- 
ever, Refs. 119). We seek to remedy this by first presenting in 
Fig. |5]the dependence of the critical current on the misorienta- 
tion angle. The three junction widths d/£, = {1-0, 1.25, 1.5} 
are chosen somewhat arbitrarily from the available range of 
values, but illustrate adequately the conditions for appearance 
of O-vr-oscillations. First, we confirm the obvious fact that the 
critical current should be symmetric with respect to (p — tt. 
Therefore we will from now on only consider the interval 
(/> e [0, tt], the maximum of the misorientation angle being the 
antiparallel orientation. Next, we see that no O-vr oscillations 
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FIG. 3: (Color online) Plot of the width-dependence of the critical 
current for several values of the spin-flip scattering rate, which is here 
taken to be uniaxial in spin space. We have defined the dimensionless 
parameter g = l/(rsfAo) as a measure of the spin-flip scattering 
rate. At d/^ — 1.0, the system is in the state. 



appear upon varying the misorientation angle for d/^ = 1.0 
fixed, which agrees with the observation from the previous 
figures that the junction appears to be in the state for all an- 
gles at this junction width. For d/£^ = {1.25, 1.5} however, 
the junction starts out in the tt state for the parallel orientation, 
and we can see that a transition takes place to the state for 
some angle ip < tt/2. This is in agreement with the result of 
Ref. [19 that a nonmonotonic dependence of the critical cur- 
rent on (p occurs when the tt state is the equilibrium state of 
the junction for = 0, and a similar statement was also made 
in Ref. Tl. 

A small effect of spin-flip scattering which may be men- 
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FIG. 4: (Color online) Plot of the width-dependence of the critical 
current for several values of the misorientation angle with spin- 
flip scattering turned off for simplicity {g — 0). 

tioned here, is that for > it may give the appearance of 
an effectively lowered misorientation angle with regard to the 
shift in the d values for O-tt crossover. As a result, for angles 
just above the critical angle, an increase of g may trigger the 
transition from complete absence of the tt state to O-tt oscilla- 
tions. 

The evolution of the critical current for variable d and cj) 
described in the previous paragraphs can be condensed to the 
phase diagram shown in Fig. |6] If only the sign of the Joseph- 
son current is of interest, each of the plots of Ic{d) may be 
thought of as a horizontal sweep through the phase diagram 
for some fixed (/>, while every plot of Ic (</>) is represented by a 
vertical sweep for some fixed d. As seen, the 7r-state becomes 
impossible above a critical angle (j>c for the present case of 
T ^ 0. We will contrast this with the finite-temperature case 
in the next section. 



B. Finite temperature 

We proceed by considering the dependence of the Joseph- 
son current on the temperature. The superconducting elec- 
trodes were assumed to be conventional superconductors un- 
affected by the ferromagnetic layers, so the standard BCS 
temperature dependence of the superconducting gap will be 
employed: 

A(r) = Ao tanh (l.74:^/TJT~^^ . (26) 

To illustrate how the critical angle 0c for 0-7r-oscillations 
depends upon increasing the temperature, including also the 
dependence on the exchange field, we plot in Fig. |8]the crit- 
ical angle as a function of the exchange splitting h for sev- 
eral values of T. As seen, 0c remains less than tt/2 up to 
H/Aq ~ 100 in the T ^ limit. However, increasing 
the temperature only slightly to T/Tc = 0.2, we see that 
(f>c rapidly approaches 7r/2. The trend is the same upon in- 
creasing the temperature even further, indicating the limit of 




FIG. 5: (Color online) Plot of how the critical current is affected 
by a change in the relative orientation of the magnetizations for 
several values of the spin-flip scattering rate, which is here taken to 
be uniaxial in spin space. 



(j)c — 7t/2 for arbitrarily high values of /i as T ~ Tc. We con- 
clude therefore that the critical angle above which O-tt oscilla- 
tions cease to exist is equal to 7r/2 as long as the temperature 
is high {T/Tc ~ 1). However, for low temperatures and weak 
exchange fields, we find that 0c deviates noticeably from 7r/2. 

Also by varying the temperature parameter may O-tt oscil- 
lations be found in the system, as shown in Figs. Q This fol- 
lows as a natural result if the critical values of junction width 
at which O-tt transitions were found in the preceding section 
are temperature-dependent. One difference from the plots of 
Ic as a function of d is the existence of no more than one tran- 
sition point (for each value of g), before the Josephson current 
necessarily vanishes at T = Tc. For increasing 0, this transi- 
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FIG. 6: Phase diagram of the Josephson junction for zero temper- 
ature, showing the regions occupied by the state and the tt state 
in width-misorientation parameter space. For the region given by 
4> £ [7i"/2, 37r/2], the n state is completely absent. We have set the 
spin-flip scattering strength to zero, g = 0. 



tion point moves leftwards until vanishing at T = at some 
critical angle. By considering the dependence on misorien- 
tation angle more carefully, we found that this critical angle, 
over which the O-tt oscillations disappear for 1^ as a function 
of T, differs from the corresponding critical angle for Ic as 
a function of d. This may be explained by going back to the 
phase diagram in Fig. |6]and noticing that (for T/Tc ~ 0) the 
maximum (j> for the tt phase region corresponds to a junction 
width d/^ « 1.35 (up to periodicity). For any other junction 
width, e.g. rf/^ = 1.25 as used in the figures, the tt state will 
be unrealizable at T = for even lower values of (f>. This 
results in the inequality that 0c for thermally induced oscil- 
lations is always less than or equal to (pc for width-induced 
oscillations. 

Another point on which the thermally induced O-tt oscil- 
lations differs from those obtained by varying the junction 
width or the misorientation angle, is the remarkably stronger 
dependence on the spin-flip scattering rate. Increasing g shifts 
the transition point significantly to the left, and furthermore 
strongly influence the ratios /c(g > 0)//c(0). Similar find- 
ings were reported in Ref. 47. 

So far, we have not considered the dependence of the crit- 
ical current on the misorientation angle while simultaneously 
going away from the limiting case of T = 0. In principle 
the phase diagrams presented might readily be generalized 
to a 3-dimensional d-<j)-T phase diagram, but we justify the 
omittance of this by arguing that the phase diagram of the 
system does not contain many interesting new features not 
already contained in the 2-dimensional projection presented 
here. However, as is clearly visible in the phase diagram for 
T/Tc — 0.5 as shown in Fig. |9] the critical angle is exactly 
0c = 7i'/2, in full agreement with the analysis done in Ref. 1. 
Another development as T increases is the 7r-region obtaining 
a more symmetric shape, also this in agreement with Ref. Jj. 



C. Enhancement effect 

As was seen in Fig. |5] there is a significant difference be- 
tween the current in the parallel configuration = and the 
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FIG. 7: (Color online) Plot of the temperature dependence of the crit- 
ical current for several values of the spin-flip scattering rate, which is 
here taken to be uniaxial in spin space. The junction width is given 
by the value d/^ = 1.25. 



antiparallel = tt. The ratio between these two critical cur- 
rents, /c(0 — 7r)//c(0 = 0), is plotted as a function of the 
junction width in Fig. [TO] The fact that one always observes 
O-TT oscillations for varying d in the parallel case, but never 
in the antiparallel, leads to a divergence of the ratio at certain 
values of the junction width, since /c(0) drops to zero at this 
transition point while Ic (tt) remains finite for all values of d. 
We note also that one will always have /c(7r) > /c(0), but the 
critical current is a monotonously increasing function of up 
to = TT only in the case that the system is in the state for 

0=0. 

This enhancement of the Josephson current by switching 
the direction of magnetization may possibly be utilized in a 
device for controlling the magnitude of the current, if the junc- 
tion is tuned to the vicinity of a transition point. A similar 
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FIG. 8: (Color online) Plot of the critical angle at which the 0- 
TT oscillations disappear, as a function of the exhange splitting h, 
for a number of temperatures. Spin-flip scattering is neglected for 
simplicity, g = 0. 



- 0.25 




FIG. 9: Phase diagram of the Josephson junction for temperature 
T/Tc — 0.5, showing the occupation of the state and the n state in 
width-misorientation parameter space in a similar manner as Fig. [5] 
does for zero temperature. For the region given by G [7r/2, 37r/2], 
the n state is completely absent. We have here set the spin-flip scat- 
tering strength to zero, g = 0. 




FIG. 10: (Color online) Plot of the ratio between the critical current 
in the antiparallel [Ic{(j} ~ tt)] and parallel [/c(<?f> = 0)] orientation 
as a function of width d. We have set g — 0, but the behaviour is 
qualitatively identical for g > 0. 



effect was mentioned by Golubov et al}'' , who considered the 
exchange field h as the variable parameter, but to our knowl- 
edge it has not yet been pointed out how this effect may be 
applied by tuning the junction width or the temperature. 



IV. DISCUSSION 

Above, we have neglected the spatial depletion of the su- 
perconducting order parameter near the S/F interfaces. This 
approximation is expected to be excellent in the case of a 
low-transparency interface."*^ Moreover, it is well-known that 
a magnetic flux threading in a Josephson junction in general 
gives rise to a Fraunhofer modulation of the current as a func- 
tion of the flux4^ We here neglect this modification by as- 
suming that the flux constituted by the ferromagnetic region 
is sufficiently weak compared to the elementary flux quan- 
tum. This is the case for either a small enough surface area or 
weak enough magnetization, but neither of these preclude the 
possibility of having an appreciable energy exchange splitting 
between the majority and minority spin bands. 

In the limit of antiparallel orientation, the tt state will be- 
come disallowed because the effect of the ferromagnetic lay- 
ers cancels, effectively giving a S/N/S junction. However, re- 
membering the symmetry requirements around = tt and the 
possibility that also partial cancellation is sufficient to ren- 
der the sign change in // o impossible, we realize that O-tt- 
oscillations may vanish for two intermediate angles (f) = ztipc 
with < 0c < TT. The partial cancellation of the exchange 
fields commences at (p > 7r/4, which may provide a clue as 
to why the critical angle is always in the vicinity of = 7r/4 
as seen in Fig. [8] Note that although the O-tt oscillations van- 
ish above the critical angle (pc, it is evident from e.g. Fig. [3] 
that the critical current does not decay monotonously even for 
(j> > (j>c as it would have in an S/N/S junction. 

One interesting observation in our study is that although the 
critical angle varies, we never find (pc > 7r/2. This means that 
any choice of parameters that brings us away from the limit 
considered in Ref. seems to lower <j)c but never increase 
it. A conjecture which may shed some more light on this 
phenomenon, is that the magnetization component of the F2 
layer perpendicular to the magnetization of the Fi layer can 
be viewed as an additional effective spin-flip scattering effect. 
It is known^^ that sufficiently strong spin-flip scattering may 
remove the oscillations in the anomalous Green's function en- 
tirely, thereby inhibiting O-tt oscillations. This effect will be 
at its strongest for ip 7r/2, which may account for the some- 
what surprising fact that (j)^ always remains close to cf) — tt/2. 
Combined with the effective cancellation of the magnetiza- 
tion described in the previous paragraph, this also serves as a 
possible explanation why the 7r-state is forbidden for the ori- 
entations tt/2 < (p < TT. To gain further understanding of 
the phase diagram of multilayer SFS junctions, we suggest to 
extend the study to a trilayer model similar as that studied in 
Ref. ^29, but where the three layers have equal thickness. If 
one fixes the middle layer and varies the orientation of the left- 
most and rightmost layers, the possible existence and value of 
a critical misorientation angle would give some hints to the 
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origin of this phenomenon also in our bilayer system. 

Finally, we would like to present a decomposition of the 
current to serve as the basis for a discussion on the long range 
contributions to the Josephson current. It can be shown easily 
from the formula for the Josephson current in Eq. ( |22] | that 
one may rewrite A/+ + = Mt — Ms where Afs and Mt 
is expressed exclusively by the components fs and /( of the 
Green's function, respectively. It can furthermore be shown 
that for the ESP components one gets = A/| = Mcr- One 
may in this manner decompose the current as 

/oo 
deRe{[-M,(e)] 

+ [Mt{e)] + [2M,(e)]} x tanh(/3e/2). (27) 

While the total current is easily shown in our framework 
to be constant throughout the junction, the separate compo- 
nents given above need not be, and the spatial dependence of 
each contribution is plotted in Fig. [TT]for selected parameter 
values. For parallel orientation, one naturally finds that the 
ESP correlations do not contribute to the current at all. For 
(j) = 7r/2, where the ESP contribution is naively expected to 
be at its maximum, we find that the OSP triplet contribution 
however, is exactly zero throughout F2. That Mt and 2AIa 
seem to change roles at a; = can be explained in a natural 
way by remembering that ultimately, the quantization axis was 
taken to be || z for all x. For a quantization axis || y, however, 
the components considered as ESP in the former case would 
here correspond to OSP components, having spin projection 

Sy = 0. 

The above argument may be used to to clarify a point re- 
garding the contribution to the current from LRTC. Ref. 5^ 
claims that a long range component of the critical current 
does not appear in a SF1F2S structure even with noncoUinear 
magnetizations, which seems at odds with our observation in 
Fig. [TT]of a non-zero ESP component to the current in F2. It 
is however important to maintain the distinction between the 
ESP contribution to the current and a LRTC contribution. As 
explained in the case of = 7r/2, the total triplet contribution 
is equal for the two magnetic layers, but appears as an ESP 
component in F2 only because of the choice of quantization 
axis. Therefore, the ESP contribution in F2 is equivalent to 
the OSP contribution in Fi, and thus cannot be regarded as a 
true long range component. A long range Josephson effect is 
defined by the absence of the exchange field in the exponent 
for the relevant Green's functions, making its decay length in 
a ferromagnet comparable to that of a normal metal. Inspect- 
ing Eqs. ( fTSb and we see that this certainly is not the case 
for = 7r/2. 

The discussion concerning the different contributions to 
the Josephson current may also hold an important clue con- 
cerning the O-TT oscillations. When inspecting the symme- 
try components separately, we observed that there can be 0- 
TT-oscillations of both the singlet and the triplet components 
simultaneously. If this is generally the case, one idea is to 
investigate for which parameters these two contributions act 
constructively and for which they act destructively. In this 
way the relative interplay of the different symmetry compo- 
nents of Ic may offer an explanation of the behaviour of the 
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FIG. 11: (Color online) Spatial dependence through the ferromagnet 
for the separate contributions to the critical current, as decomposed 
in Eq. J27t . Ic,{x) refer to Ic.s, Ic,t and Ic,<y, with contributions to 
the integrand from the components [~AIa{e)], Mt{e) and 2Ma-{e), 
respectively. The current was evaluated for a junction width given 
by d/^ — 1.0 and temperature given by T/Tc = 0.4, and spin-flip 
scattering is neglected for simplicity (g = 0). 

critical angle 0c for which the O-tt oscillations of the total crit- 
ical current disappear. 

V. SUMMARY 

In conclusion, we have investigated the O-tt oscillations 
of the critical current in a diffusive SF1F2S Josephson junc- 
tion with noncollinear magnetization, where the effects of 
noncollinearity and spin-flip scattering have been studied in 
particular The introduction of spin-flip scattering does not 
change the Josephson current dramatically, so the phase di- 
agrams presented above for zero spin-flip scattering would 
therefore not be qualitatively changed much by setting g > 0. 
Also, comparing isotropic and uniaxial spin-flip scattering, we 
found that the effect of these was very similar both qualita- 
tively and quantitatively. 

Oscillations of the critical current for varying junction 
width disappear when the relative angle of magnetization 
passes a critical value (pc, making the tt state unattainable re- 
gardless of choice of the other parameters. This critical an- 
gle equals 7r/2 in the limit of relatively high temperatures 
(T/Tc ~ 1), but is lowered below tt/2 when the temperature 
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is low and the exchange field is small simultaneously. These 
dependencies on the various parameters can rather easily be 
read out of phase diagrams of the kind we have presented. A 
straight-forward analytical approach to the behaviour of the 
critical angle is challenging due to the many variables in our 
system, but we have discussed several ways by which its ori- 
gin can be better understood. 

With regard to the effect of finite interface transparencies, 
we have not found any signs throughout our investigations im- 
plying that the results of Ref. 1 change in any significant way. 
Our mapping of the relevant parameter regimes does however 
serve as a starting point for looking for new interesting effects 
that may appear upon varying the parameters kept fixed in our 
case, i.e. the transparency 7 and the exchange field h in par- 
ticular A natural course for a continuation of this work would 
be expanding the system from a bilayer ferromagnet to a tri- 



layer, similar to the Josephson junction considered in Ref. 
where the relevant parameter regime for a significant LRTC 
contribution to the Josephson current was found. It might be 
interesting to see how LRTC manifests in our framework of 
symmetry components to Ic, and to investigate for what re- 
gion in parameter space the vr state can be realized in such a 
system. 
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